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We present recent results of B and charm decays from the BABAR experiment. These results 
include searches for rare or forbidden charm decays, measurements of | V u t>\ from inclusive B X u lv 
t-h decays , observation of the semileptonic B — > D^t~v t decays, direct CP violation asymmetry in 

B — > X s+ d^ and in D + — > K^ir^, and T-violation in Df-. — > K^Kgir + ir~ . These studies are based 
on the final dataset collected by BABAR at the PEP-II B factory at SLAC in the period 1999-2008. 
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1 Introduction 



We present some recent results concerning B and charm decays from the BABAR experiment. These 
results include searches for rare or forbidden charm decays, measurements of \V u b\ from inclusive 
B — > X U W decays, observation of the semileptonic B — > D^t~v t decays, direct CP violation 
asymmetry in B — > X s+ d7 and in D + — > K^tt + , and T- violation in D~^\ — > K + K®it + it~ . The 
aim of these measurements is to test the Standard Model (SM) mechanism of CP violation and to 
search for possible contributions from New Physics (NP) beyond the SM. 

Important areas of search for NP are processes which are expected at low level in SM and 
which could be enhanced by NP. This is the case for example of the measurement of branching 
fraction of many rare decays , in particular Flavor Changing Neutral Currents (FCNC), and the 
measurement of CP violating asymmetries. Such measurements are sensitive to NP scenarios and 
have been already very powerful in constraining the parameter space of NP models. 

The studies presented in this talk are based on the final (or almost final) dataset collected by 
BABAR in the period 1999-2008 at the PEP-II B factory at SLAC. 



2 D° ->> 77 and D° A 

FCNC decays, forbidden at tree level in SM [I] but allowed at higher order, have been already 
observed in K and B meson systems [2j. In the charm sector a FCNC process has an additional 
suppression due to the GIM mechanism JIJ. Thus far no charm FCNC decay process has been 
observed. Interest in FCNC processes in the charm sector increased with the recent measurements 
of D° - D° mixing [3]. Source of this mixing in fact can come from NP, enhancing the branching 
fraction of FCNC with respect to SM calculation (4j[5]. 

Calculations in the framework of vector meson dominance E] and of heavy quark effective theory 
combined with chiral perturbation theory [ij] predict for the D° — > 77 decay the dominance of long- 
range effects and a branching fraction of (3.5^2 g) x 10 -8 and (1.0±0.5) x 10~ 8 , respectively. These 
branching fraction estimates are orders of magnitude below the sensitivity of current experiments. 
But gluino exchange in Minimal Supersymmetric Standard Model (MSSM) can enhance the SM 
branching fraction up to a factor 200 [5] (so within BABAR sensitivity). 

CLEO has measured an Upper Limit (UL) of 2.9 x 10~ 5 at 90% Confidence Level (CL) for the 
branching fraction of the decay D° — > 77 |7| and a branching fraction of (8.1 ± 0.5) x 10 -4 for the 
decay mode D° — > -k°tt° (§]. 

Recently BABAR has studied the decays Z) u — > 77 and D° -» vr 7r° with a dataset corresponding 
to an integrated luminosity of 470.5 fb^ 1 . D° is reconstructed with a D* + tag in the decay 
D* + — > Z)°7r + 19], with the charge of the soft pion from D* + indicating the initial flavor of the D°. 
This tag suppresses the dominant combinatoric background. To avoid uncertainties in the number 
of D* + , the branching fractions of the decay modes D° — > 77 and D — > 7r°vr are measured relative 
to the reference decay mode D° — > K^tt . This mode has a large and precisely measured branching 
fraction (1.22 ± 0.05) x 10" 2 [loj. 

D° candidates are reconstructed in the decays D° — > 77, D° — > 7r°7r°, and D° — > Kgir . In 
the decay D° — > 77 the photon candidates have a center-of-mass (CM) energy between 0.74 and 
4 GeV and the main background D — > TT°TT° is suppressed with a tt veto. B background in 
D° -> 77 ( D° — > 7r°7r°) is rejected selecting D* candidate with a CM momentum above 2.85 
(2.4) GeV/c. Signal (reference mode) selection efficiency for the D° — > 77 mode is 6.1 (7.6) %. For 
the D° —7- 7r°7r° mode the signal (reference mode) selection efficiency is 15.2 (12.0) %. Signal yields 
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are extracted with unbinned Maximum Likelihood (ML ) fit to the D° invariant mass distribution 
(Fig. [I]). 




Figure 1: Fit results to D° — > 77 (on the left) showing data (points), the combinatoric background 
(long-dashed red), the combinatoric background plus the D° — > 77 background (short-dashed 
magenta) and and the total fit function including signal and the two backgrounds (solid blue); fit 
results to D° — > 7r°7r° (on the right) showing data (points), the combinatoric background (dashed 
red) and the full fit function (solid blue) including signal and combinatoric background. 



For D — > 77 the fit signal yield is —6 ± 15. Including statistical and systematic uncertainties 
the UL at 90% CL for this decay mode is 2.4 x 10~ 6 . This preliminary result is in agreement with 
SM expectation. Based on this result enhancement due to gluino exchange cannot exceed a factor 
70. 

The preliminary result for the branching fraction of the D° — > 7r°7r° decay is (8.4 ± 0.1 ± 0.4 ± 
0.3) x 10 -4 , where the uncertainties refer to statistical, systematic, and reference mode branching 
fractions uncertainties, respectively. 



2.1 Search for X c ->• hl + l 



In a recent paper [ll] BABAR searched for charm hadron decays of the type — > h^Pl^ + [9], 
where X^~ is a charm hadron (D + , Df or A+), /*• ) + is an electron or a muon, and can be a pion 
or a kaon (a proton in A+ decay modes). The analysis is based on a data sample of 384 fb^ 1 of 
e + e~ annihilation data collected at or close to the T(4S') resonance. 

In total 35 decay modes have been studied. Among them there are decay modes with oppositely 
charged leptons of the same flavor which proceed through FCXC. These decays are very rare. 
There are decays with oppositely charged leptons of different flavor. These correspond to lepton- 
flavor violation decay modes which in SM are essentially forbidden because they can proceed only 
through lepton mixing. There are decays with two leptons of the same charge. These are lepton- 
number violating processes which are forbidden in the SM. The most stringent existing ULs for 
the branching fractions of the decays X^~ — > h^r^l^ + are in the range [1 — 700] x 10 -6 [12 - 16 
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Figure 2: Standard Model short-distance contributions to the c — > ul + l~ transition. 



Transitions c —¥ ul + l proceed in SM through diagrams shown in Fig. [2] and are expected with 
branching fractions of 0(1O -8 ) 4p!7 . Several extensions of the SM predict an enhancement of these 



branching fractions |4| |17[|18| . The target decay modes may have also long-distance contributions 
from leptonic decays of intermediate resonances like Du\ — > X U V with V — > l~l + which are 
expected with a branching fraction of 0(1O -6 ) [4|fl7] . At current experimental sensitivity these 
long-distance contributions can come only from -4 7r + (p decays with cj) — > The effect of 

these long-distance contributions is suppressed with a cut on the invariant mass around the 
(p meson. 

In the signal reconstruction three tracks, one identified as n, K, or proton and two identified 
as electron or muon )), are merged. Charm hadron is selected with a momentum p* in e + e~ 
CM > 2.5 GeV/c to suppress charm hadron production from B decays. The final selection is done 
with a likelihood ratio (LR) using three discriminating variables (p* , total reconstructed energy in 
the event, and flight length significance of the charm hadron candidate). The minimum LR value 
is chosen independently for each decay mode. 

Signal yields are extracted with extended unbinned ML fit to the invariant mass distributions 
of Ml. There are three components in the fit: signal, combinatoric background, and background 
from nonleptonic charm decays in which two hadrons are misidentified as leptons. The fitted signal 
yields are translated into branching fractions by normalizing them to the yields of known charm 
decays with similar kinematics (Dt^ -4 ir + (p (<ft -4 KK) and A+ -4 pK~ 7r + ). No signals are 
found and Bayesian ULs at 90% CL are calculated. Limits in 32 of the 35 studied decay modes are 
improved upon the previous ones, in most cases by more than order of magnitude. 

3 Measurement of \V u h\ from Inclusive B — > X u lv Decays 

The magnitude of the Cabibbo-Kobayashi-Maskawa quark-mixing matrix element V u b UM can be 
determined from inclusive semileptonic B decays to charmless final states X U W, where I = e or a, 
and X u is a hadronic system (without charm). The real difficulty in this inclusive measurement 
comes from the overwhelming charm background from B — > X c lv which has a rate fifty times larger 
and an event topology very similar to signal. 



In a recent analysis 20 BABAR, using the full dataset of 467 x 10 6 BB pairs, has measured the 
partial branching fractions of the B — > X u lv decays [9] , restricting the analysis in selected regions 
of the phase space where most effective is the suppression of the charm background. The event 
selection uses a hadronic tag: in the sample of T(45) — > BB events, one B decaying into hadrons 
is fully reconstructed (B reco ) and the other B (B recoi i) is identified by the presence of an electron 
or a muon. The B reco is reconstructed in many exclusive hadronic decays B reco — > D^Y ± , where 
D^*' is a charmed meson (D°, D + , D*°, or D*^ ) while the charged hadronic system consists 
of up to five charged hadrons, pions or kaons plus up to two neutral mesons (Kg, ir°). 
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In the B Teco n rest-frame we require only one charged lepton with momentum > 1 GeV/ c and 
the hadronic system X u is reconstructed from charged particles and neutral clusters not associated 
to the B reco or to the charged lepton. Neutrino is reconstructed from the missing four-momentum 
in the whole event. Requirements on several kinematic observables were applied in different phase 
space regions to select the final signal events. 

Partial branching fractions are measured in several regions of phase space and are normalized 
to the total semileptonic branching fractions, thus reducing several systematic uncertainties. Con- 
sidering the most inclusive measurement (based only on the requirement p* > 1.0 GeV/c), from a 
two-dimensional Mx — q 2 ML fit to the hadronic invariant mass Mx and the leptonic mass squared 
q 2 , we measure: 

AB{B ->■ X u lu;p* > l.OGeV/c) = (1.80 ± 0.13 ± 0.15 ± 0.02) x 10~ 3 , (1) 

where the first uncertainty is statistical, the second systematic, and the third theoretical. 
The measured partial branching fractions are related to \V u b\ via the following relation: 



I Kb I 



>AB(B -»• XJv) 

TflAr \heory 



where tb is the B meson lifetime and AT t heory is the theoretically predicted AB(B — > X u lv) 
partial branching fraction for the selected phase space region. This prediction is calculated on 
the basis of four different QCD models: Bosch, Lange, Neubert, and Paz (BLNP) |22| , Gambino, 
Giordano, Ossola, and Uraltsev (GGOU) [23], Andersen and Gardi (DGE) [24], and Aglietti, Di 
Lodovico, Ferrera, and Ricciardi (ADFR) [25] . Results for \V u b\ for these four different QCD 
calculations using the partial branching fraction obtained in the most inclusive measurement (Eq. [T]) 
are presented in Table [T] In the last row of Table [T] we give the arithmetic average of the values 
and uncertainties obtained with the four QCD calculations (the first uncertainty is experimental 
and the second theoretical). The total uncertainty is about 6.9 %, comparable in precision with 



the Belle result 26 



Table 1: Results for \V u b\ for the four different QCD calculations for the most inclusive partial 
branching fraction measurement. 



QCD Calculation 


\v ub \(w- 3 ) 


BLNP 
DGE 


4.27 ± 0.15 ±0.18+0:20 
4.34 ± 0.16 ± 0.18±g;f| 


GGOU 
ADFR 


4.29 ± 0.15 ±0.18l{j:i4 
4.35 ±0.19 ±0.20+^1 


Arithmetic Average 


4.31 ±0.25 ±0.16 



The value of \V u b\ obtained in this inclusive analysis is higher than the value obtained by BABAR 
in an exclusive analysis [21]. The discrepancy is at a level of about 2.7 a and is increased in the 
latest measurements obtained with significantly decreased uncertainties thanks to the improved 
experimental techniques and theoretical inputs. A similar discrepancy is also present in BELLE 
results [26l[27 . 
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4 Measurement of the decays B — > Dt v t and B — > D*t v t 



28 



30 



Semileptonic B decays to r lepton B — > D^t~u t |9| can provide constraints on the SM 
and are sensitive to physics beyond the SM. In extensions of the SM (such as in the multi- 
Higgs doublet models and the MSSM), intermediate charged Higgs boson can contribute to the 



amplitude, modifying significantly the branching fraction 31-36 . Branching fractions of these 
decay modes are smaller compared to those of decay modes to final states containing a light lepton, 
1 = e or fj,. SM expectations for the relative rates between signal and reference modes are: 7£(D) 
= B(B -> Dtu t )/B(B -»• Dlvi) = 0.31 ± 0.02 (3T] and K(D*) = B(B -»• D*tv t )/B(B -> D*lv{) = 
0.25 ± 0.02 [35]. Multi-Higgs doublet models predict an effect on the ratio 7£(D) much stronger 



than on 1Z(D*) 31-36 . This effect may enhance or decrease TZ(D^), depending on the value of 



the ratio tan /3 / m^± of the Higgs parameter tan (3 and the charged Higgs mass m H ± . 
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Recently 38 BABAR updated with the full data sample of 471 x 10 BB pairs previous analyses 
of the decays B — > D^*'tv t . T(AS) — > BB events are tagged by the hadronic decay of one of 
the B mesons (B tag )- The tag decay modes are of the type B tag — > SX^, where S can be D, D*, 
D s , D*, or J/ip reconstructed in many s. is a charged state with a maximum of 5 particles (tt 
or K) including up to two neutral particles (tt° or Kg). 

For each Bt ag candidate in a selected event, the other B{B s i g ) is searched for combining a single 
charged lepton and a meson. The r lepton is reconstructed only in the purely leptonic decays 
r~ — > e~v e v T and r~ — > h^v^Ut while the £)W meson is reconstructed in the four modes D°, D*°, 
D + , and D* + . D° (D + ) is reconstructed in 5 (6) decay modes for a combined branching fraction 
of 35.8% (27.3%). D* meson is identified in the decays D* + D°tt + , D + tt° and D*° -> D°ir°, 
D°j. The signal modes have in the final state one secondary lepton and three neutrinos while the 
reference modes have a primary lepton and one neutrino. 

For the separation of signal and reference modes the most discriminating variable is the missing 
mass squared , defined as vr? mif , s = (p e + e - — Ptag — Pd(*) ~ Pi) 2 > where p are four-momenta. This 
quantity peaks at zero in the reference decay modes where only one neutrino is missed while in the 
signal modes the vr? miss distribution is broad and extends up to 8 (GeV/c 2 ) 2 . Further separation 
of signal and reference modes is obtained with a requirement on the minimum momentum transfer 
. For decays with a t, q 2 = (p T + p UT ) 2 > m 2 . = 3.16 (GeV/c 2 ) 2 . The applied requirement q 2 > 4 
( GeV/c 2 ) 2 retains 98% of the signal decays and rejects more than 30% of the reference modes. 



Table 2: Results for the ratios TZ(D^) , the individual signal branching fraction, and the sig- 
nal significance Stot including systematic uncertainties, and the significance S s tat (only statistical 
uncertainties) 



Decay Mode 




B(B -> D^ru) (%) 


Stot (Sstat) 


D u t~V t 


0.422 ±0.074 ±0.059 


0.96 ±0.17 ±0.14 


5.0(6.2) 


D*°t-u t 


0.314 ±0.030 ±0.028 


1.73 ±0.17 ±0.18 


8.9(11.9) 


D+t-V t 


0.513 ±0.081 ±0.067 


1.08 ±0.19 ±0.15 


6.0(7.5) 


D*+t-V t 


0.356 ±0.038 ±0.032 


1.82 ±0.19 ±0.17 


9.5(12.1) 


Dt~V t 


0.456 ± 0.053 ± 0.056 


1.04 ±0.12 ±0.14 


6.9(9.6) 


D*t~V t 


0.325 ±0.023 ±0.027 


1.79 ±0.13 ±0.17 


11.3(17.1) 



Combinatorial background is reduced constraining to the same vertex the charged daughters 



6 



of the flW and B mesons. Improved discrimination of signal and reference modes from several 
backgrounds is obtained using a multivariate discrimination with a Boosted Decision Tree (BUT): 
12 discriminating classifiers using 8 discriminating variables each. 

Source of difficult background is the decay B — > D**l~i>i. D** mainly decays to D^ir which 
enter in the selection when the pion is neutral and not reconstructed or is charged and associated to 
Btag- This background has been studied with four control samples (D(*'tt°) identical to the signal, 
except for an additional ir° selected in the mass range [120 , 150] MeV. 

To extract yields for signal and reference modes we perform an extended, unbinned ML 2D 
fits to the distribution of vni^ niss vs \p\\ (lepton momentum calculated in the B meson rest frame). 
Simultaneous fit to the four signal channels (D°, D*°, D + , D* + ) and the four D^ir channels. 
Another fit is performed imposing the isospin constraint: K(D°) = K(D + ) = TZ(D) and K(D*°) = 
K{D*+) =K(D*). 

The fit describes reasonable well the four D^*'w° channels inside the sizable statistical uncer- 
tainties. The fit describes very well the large contributions of the reference decay modes. Fit 
results for D*° and D* + channels are very good. Both channels are observed with a significance 
higher than lie (only statistical uncertainties). For the D channels the fit projection onto m'^ niss 
shows an excess of data in the range 1.5 < Tri^ss < 3.5 (GeV/c 2 ) 2 and an overestimate of events 
for m^ niss > 5 (GeV/c 2 ) 2 . These regions are dominated by continuum and B combinatorial back- 
grounds which are fixed in the fit to what is expected by simulation. It is not clear if the fit-data 
differences are statistical or systematic. Both D° and D + channels are observed with a significance 
(statistical uncertainties only) higher than 6 a. 

Main sources of systematic uncertainties are the MC simulation of the background, the statis- 
tical uncertainties of the simulated samples, and the D**lv decay modes. In order to understand 
the origin of the difference in the fit-data comparison, the BUT requirements have been changed 
in such a way to have a sample passing the selection 50%. 80%, 120%, and 200% compared to 
the nominal sample. The agreement between fit and data improved using both more and less re- 
strictive BUT requirements. We assign a systematic uncertainty equal to half of the variation on 
TZ(D^ in the fit when applying tight BDT requirements (50% nominal sample) and loose BUT 
requirements (200% nominal sample). This systematic uncertainty (9.5% on TZ(D) and 6.5% on 
1Z(D*)) is comparable in size with statistical uncertainties. It should be eliminate or reduced once 
the source of the difference in fit-data is understood. 

Table [2] summarizes the results obtained from the two fits: the first in which all four signal 
yields can vary independently, and the second (last two rows in the table) in which isospin relations 
are imposed. 

These preliminary results are in agreement with previous BaBar results [39] and with Belle 
measurements 40 and have significantly reduced uncertainties. The decays B° — > D + t~v t and 



B —7- D°t v T are observed for the first time. 
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Figure |3[shows 11(D) as a function of the ratio tan /3 / mjj± . The violet band is the theoretical 
prediction |41| while the horizontal blue band is the present measurement. Present result favors 
large values of tan/3/m#±. Furthermore TZ(D) is about 1.8 a in excess over SM prediction. A 
similar excess is also measured for 1Z(D*). 
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Figure 3: 71(D) as a function of the ratio t&n(3/m H ±. 



5 Direct CP Asymmetry in Inclusive B — > X s+( pj 

In the SM the inclusive electromagnetic radiative decays b — > s*y or b — > proceed at the leading 
order Feynman diagram shown in Fig. |4j 




U 5 C,t 

Figure 4: Leading order diagram for b — > (s,d)^. 



The branching fraction of the process B(B -)■ X s -f) = (3.15 ± 0.23) x 10" 4 with £ 7 > 1.6 GeV 
(in the B meson rest frame) has been calculate in SM at next-to-next-leading order with a precision 
of 7% [?]. Because new heavy particles may enter in the loop at leading order, the value of this 
branching fraction is highly sensitive to NP 43 . Another quantity highly sensitive to contributions 



of NP is the direct CP asymmetry [44] defined as: 



A C p 



T(b -> S7 + b -> dry) - T(b ->• s-f + b -)■ dq) 



T(b-> sj + b^ dj) + r 
This asymmetry in SM is expected ~ 10 -6 



S7 + b — > c?7) 



44 



with the asymmetries for b — > sj and 6—7-^7 



opposite with nearly exact cancellation. In NP scenarios Acp can be at about 10% (45). 
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In a recent analysis BABAR has measured direct CP asymmetry in inclusive B -4 ^ s +d7 using 
a data sample of 383 x 10 6 BB pairs. The reconstruction of the event uses a fully inclusive method 
with a semileptonic tag (Fig. [5]) . 



Semileptonic 
tag 

e" or \i~ 




Figure 5: Semileptonic tag. 



Signal events are identified by high energy photon and the X s is not reconstructed. Therefore 
is not possible to distinguish X s and Xj, states and what is measured is B — > X^ s+d ^. The 
Btag is searched for in a semileptonic decay mode, checking that the remaining particles in the 
event are consistent with a B decay. Photon selection requires at least one photon with energy 
1.53 < E* < 3.5 GeV and the tag lepton can be an electron or a muon with momentum p* > 1.05 
GeV (E* and p* are calculated in the T(45) rest frame). Continuum background is suppressed 
using a neural network discriminant based on eight topological variables. Remaining continuum 
background in the final sample is estimated using off-resonance data collected 40 MeV below the 
T(4S) resonance. The BB background is mostly due to photons from low- mass mesons (mainly 
7T° and r]). This background is removed using explicit vetoes. Remaining BB background is 
estimated from MC simulation and is cross-checked against data with control samples. Continuum 
and BB background estimation have been validated with control samples provided by the photon 
spectrum. 

In Fig. [6] we show the signal and backgrounds distributions with only the high energy photon 
requirement (on the left in logarithmic scale) and after all selection requirements (on the right in 
linear scale). 




Figure 6: Photon spectrum after requiring a high energy photon on the left and after all selection 
requirements on the right. 
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Direct CP asymmetry is insensitive to the photon energy cut. So to reduce the sensitivity to 
background the asymmetry is calculated in the optimized range (2.1 — 2.8) GeV. Lepton charge 
gives B flavor and separation of B and B is done according to the charge of the lepton tag. The 
tagged signal yields are N + = 2623 ± 158 and N~ = 2397 ± 151 and the measured CP asymmetry 
is: 



N+ - N~ 
N+ + N- 



A^ as = — — = 0.045 ± 0.044 



This result must be corrected for the dilution due to mistag fraction u>, Acp = > f° r the 

uncertainty in the BB background estimation, and for the bias induced by CP asymmetry in the 
BB background and by detection asymmetry. Correcting AQp as for these effects we obtain the 
preliminary result: 

A C p = 0.056 ± 0.060 stat ± 0.018 S2/St 

No significant asymmetry is observed in agreement with SM expectations. This result is the 
most precise to date |46|. 



6 Search for CP Violation in D + -» K° s tt + 

CP violating asymmetries have been measure both in K and B systems with results in agreement 
with SM expectations. CP violation has yet not been observed in charm decays, where the SM 
expectations for CV violating asymmetries are at the level of 10 -3 or less [47 j. 



In a recent analysis 48 BABAR searched for CP violation in the decay D ± — > K^ir^ , measuring 

sr: 

r(r>+ — ± k®.tt+\ - v(r>- k^-) 

A C p 



the direct CP violating parameter: 

T(D+ -> K%n + ) - T(D~ -> K° s ir- 



F(D+ -> K®tt+) + T(D- -> Kgir~) 
with T partial decay width of the decay. 

Although the SM prediction for direct CP violation due to the interference between Cabibbo- 
allowed and doubly Cabibbo-suppressed amplitudes is negligible [49], K°- K° mixing induces a 
time-integrated CP violating asymmetry of (— 0.332±0.006) % (lOj. Contributions of physics beyond 



the SM may enhance the value of Acp up to 1% 49,501. Previous measurements of Acp have been 



reported by CLEO-c [-0.6 ± 1.0 ± 0.3 %) [5l] and Belle Collaboration [-0.71 ± 0.19 ± 0.20%] [52] 
We reconstruct Df — > K^tt^ 1 decays combining a Kg candidate with a charged pion candidate. 
The Kg decays to tt + tt~ and is selected with an invariant mass within ±10 MeV/c 2 of the nominal 
Kg mass [10 . The reconstructed pion candidate is selected with a momentum px in the plane 
perpendicular to the z axis greater than 400 GeV/c. The selected candidate has an invariant 



mass within ±65 MeV/c 2 of the nominal D + mass 10 , and a momentum p* in the e + e~ CM in 

the range (2 — 5) GeV/c. Further suppression of background has been achieved using a BDT with 

seven event topological discriminating variables. 

Signal yield is extracted with a binned ML fit to the invariant mass distribution of the selected 
candidates. There are three fit components: signal, a background from D^r — > KgK^ and a 

combinatoric background. 
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1.82 1.84 1.86 1.88 1.9 1.92 

m(K§jr*) [CeV/c 2 ] 



Figure 7: Invariant mass distribution for Kgir : solid curve is the fit to the data (points), the 
dashed line is the sum of backgrounds, and the dotted line is the combinatorial background only. 
The vertical scale on the plot is logarithmic. 



Data and fit are shown in Fig. [7| The signal yield asymmetry is: 

, = N D+ - N D - 
N D+ + N D - ' 

where N D + and N D - are the fitted yields for D + — > K^it + and D~ — > Kgir~ , respectively. 

The quantity A includes contribution not only from Acp but also from the forward-backward 
asymmetry Afb m e + e~ — > cc due to 7* — Z° interference and other QED processes. Another 
source of asymmetry (A e ) contribution to A is induced by the detector as a consequence of the 
difference in reconstruction efficiency of D + — > K^ir + and D~ — > Kgir~ due to differences in 
reconstruction and identification efficiencies for 7r + and 7r~. So we can write: 



A = A C p + A FB + A e (2) 

A e has been measured using a control sample of BB decays. The bias of +0.05% to Acp has 
been included in the systematics. 

We separate Acp and Afb in Eq. [^considering that Afb is an odd function of cos 9* D , where 9* D 
is the polar angle of the candidate momentum in the e + e~ CM, while Acp is an even function 
of cos 0* D . Therefore the two asymmetries Acp and Afb can be written as a function of | cos#£,| 
as: 

^BdcOS^I) = A( + |COS ^ l) / ( - |C ° S ^ l) 



A C p{\ cos e* D \ 



A{+\cos6* D \) + A{-\cos6* D \ 



2 

The selected sample is divided in subsamples corresponding to 5 bins of | cos6* D \ and a simul- 
taneous binned ML fit is performed to extract signal yield asymmetries. 
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The five measured values of the parameters Acp and Afb ar e shown in Fig. [8j We measure: 

A CP = (-0.39 ± 0.13 ± 0.10)% , 

where the first uncertainty is statistical and the second systematic. This value is in agreement with 
the prediction of the SM. 
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Figure 8: On the left Acp with the central value (solid line) and the ±1<t interval (hatched region), 
and Afb (on the right) in bins of | cos^j. 



7 Search for T Violation using T-odd Correlations in Df^ — > K + K s ti + ti 

T violation in the channels D + — > K + Kgir + ir~ and Df — > K + Kgir + 7r~ [9] can be measured using 
T-odd correlations. Using the momenta of the final particles in the rest frame, the T-odd 
correlation observable can be written as: 

C T =PK+ ■ (Ptt+ XPtt-) 

This triple product is odd under time reversal. Assuming CPT invariance, T violation is 
equivalent to CP violation. We measure the asymmetry: 

_ T(C T > 0) - T(C T < 0) 
T ~ T{C T > 0) + T{C T < 0) ' 

where T is the decay rate of the process. 

Because Final State Interactions (FSI) can produce an asymmetry At different from zero |53| , 
we measure also the T-odd asymmetry for the CP-conjugate decay process: 

-j = T{-C T > 0) - T{-Ct < 0) 
T = T(-C T > 0) + T(C T < 0) ' 

where Ct = Pk- ' (Pn- x Pn+) with momenta calculated in the rest frame. To remove FSI 
effects we measure the quantity At- 

At = ^(A t -A t ) (3) 
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which is an asymmetry characterizing T violation in the weak process [54] . 

In a previous BABAR analysis [55] done on the neutral decay D° — > K + K~tt + tt~ no evidence 
of CP violation has been found. 

We describe here the search for CP violation using T-odd correlations in the decays D + — > 
K + Kgir + TT~ and Df — > K + Kgir + ir + [56]. In the reconstruction we have considered inclusive D 
decays (e + e~ — > XDt s ^), selected using kinematic constraints and particle identification. In the 
vertex fit, requiring that all three tracks with net charge +1 originate from a common vertex, a 
X 2 fit probability PI > 0.1 % has been imposed. A second fit constraining the three tracks to 
originate from the e + e~ interaction region was also done. In this case the x 2 fit probability P2 is 
large for most background events with tracks originating in the interaction region. 

D^s candidate must have a momentum p* in the CM greater than 2.5 GeV/c. To improve 
signal and background separation we consider the signed transverse decay length: 

d-pT 
\Pt\ 

where d is the distance vector between IR and the Z)jK decay vertex in the transverse plane 
and pt is the transverse momentum of Dfy With the variables p*, P1-P2, and Lt we construct 
a LR to optimize the signal yields separately for D + and Df . Separation requirements where 
optimized maximizing the statistical significance S/y/S + B, where S is the number of signal and 
B the number of background events in the signal region. 

Then the dataset is divided in four samples depending of the charge and sign of Ct 

(Ct) and on these four datasets a simultaneous fit to mass spectra has been done to extract the 
asymmetry parameters At and At- The measured values are: 

A T {D + ) = (+11.2 ± 14.1 ±5.7) x 1(T 3 

At(D') = (+35.1 ± 14.3 ± 7.2) x 1(T 3 

A T (D+) = (+99.2 ± 10.7 ± 8.3) x 10~ 3 

At(D-) = (+72.1 + 10.9 + 10.7) x 10~ 3 , 

where the first uncertainty is statistical and the second systematic. FSI Effects are larger in Df 
than in D + decays. A study of such effects can be found in ref. |57| . The T violation asymmetries 
obtained using Eq. [3] are: 



A T {D + ) = (-12.0 + 10.0 + 4.6) x 10~ 3 
A T {Df) = (-13.6 + 7.7 + 3.4) x 10~ 3 

T violation parameter is consistent with for both the two decay modes. 



8 Summary and Conclusions 

We have presented BABAR results of recent searches for rare or forbidden charm decays , mea- 
surement of the magnitude of the CKM matrix element V u b, the observation of the semilep- 
tonic B — > D^rv T decays, searches for direct CP violation asymmetry in B — V and 
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in D + —7- Kgir + , and T- violation in — > K + Kg7r + ir~ . All these analyses use the final BABAR 
dataset. Results have been improved over previous analyses and are all consistent with SM expec- 
tations. No CP violation has been observed in c — > s transitions , neither from SM nor from NP. 
More stringent limits are set on space parameters of NP models. 

B-factories reached their sensitivity limit. We expect soon a significant impact in flavor physics 
from LHCb experiment |58| and (in a few years) from the next generation super B factories (SuperB 



59 at the Cabibbo Laboratory (Tor Vergata, Rome) and SuperKEKB 60 at Tsukuba. 
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